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Abstract. To our knowledge there are no complete results expressed in terms of eigenfunctions (even not strictly 
proved mathematically) related to the system of three or more charged quantum particles. For the system of the 
three such identical particles we suggest the asymptotic formula describing the behavior of eigenfunctions at infinity 
, in configuration space. 

Introduction. The principal questions of the scattering theory in the system of two, three and more quantum 
particles, considered in terms of eigenfunctions (EF), have been studied relatively completely [l but not exhaustively. 
In this field of problems there is a set of important systems which so far remain unexplored. Here we mean the 
systems of three and more particles interacting by pair coulomb potentials. There are two widely practised (for 
computer calculations) approaches to such problems. The first one consists in the regularization of the coulomb 
potential by its substitution by for example the Yukawa potential. Until now, however, there was no possibility to 
estimate the effect of such a regularization, as no reliable (even not strictly proved) results for many-body systems 
with the pair coulomb interactions were known. 

The second approach consists in the application of the so called BBK-approximation to the continuum spectrum 
EF, refer [2]. We will describe it in detail later. In a few words, BBK-approximation can be presented by the explicit 
formula which has a discrepancy in the equation, decreasing at infinity quickly enough. The problem is that on 
some specific multi-dimensional directions at infinity, which we will call the screens where the particles approach each 
other in pairs, the discrepancy loses a satisfactory pace of decrease. Our result consists in the modification of BBK- 
: approximation on the screens, so that the modified formula allows a satisfactory description of the EF asymptotics 
at all directions at infinity. Thus, we can say that we have succeeded for the first time in deriving the asymptotic 
behavior of EF. To estimate the importance of this result, note that when considering the EF even for the quickly 
decreasing potentials, we always start the constructions with the derivation of the asymptotics of these EF (or, which 
is almost equivalent) with the derivation of the corresponding integral equations. 

Weak asymptotics. We consider the system of three three-dimensional particles of equal masses, interacting 
■ by identical pair coulomb potentials. The assumption about equality of masses and potentials is introduced only for 
', the simplicity of the presentation. These restrictions can easily removed. It is quite obvious that the system has a 
pure continuous and infinitely-tuple degenerate spectrum spanning the half-axis [0,oo). The traditional way of the 
identification of the EF is connected with the description of its possible behavior at infinity in the configuration space 
of the system. For the quickly decreasing at infinity pair potentials such an identification of the EF consists for 
example in the separation (as the main order) of the plane wave with the fixed wave vector from its asymptotics. For 
the the case of the coulomb potentials the plane wave can not be considered as the main order of the asymptotics of 
the EF in conventional sense. 

Even for the case of one particle scattering on the coulomb potential v, u(x) = ^, x e R^, x — ^, such a simple 
separation of the main order term in the asymptotic behavior of the solution is impossible. The simplest description 
can be done in terms of the so called weak asymptotics [S]. It can be presented by the formula [B]: 

V'(x,k) ^ ^ (<5(x,-k)e-^^-+"'"- - s,(x,k)e^'=--''''"-) , ^ = |x| ^ oo. (1) 

Here the coefficient Sc, in fact the scattering matrix, in the problem set up process remains undefined and is determined 
only by the complete solution of the problem. The main feature of this asymptotic description is that we consider it in 
terms of distributions relatively to the x. In spite of the weakened character of the condition, it successfully separates 
the unique solution, standard EF. 
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To make it easier for the reader to compare this formula with the classical assymptotic description of the EFs, note 
that in the weak sense e*<^'''> ^ ^ (<5(x, — k)e^*'^^ — (5(x, k)e''''^^ . 

Turn to the system of three particles, fix the motion of the center of mass and will describe the internal degrees 
of freedom in the terms of the variable z (details follow below). Then the assymptotics of EF ^(z,q), q - is a wave 
vector, which can be considered as a weak alternative of the plain wave for the coulomb potentials, characterized by 
the formula 



*(z,q)^2(^— j '5(q,-z)e ^^^^'^^^ 

The system of descriptions for the vectors reproduces here the descriptions used in the formula ([T]), only with the 
difference that z G F, diniF = 6. Here again S'c - is a function indefinite when setting up a problem, which as a final 
result must turn into a scattering matrix. 

We definitely believe that such a characterization of the EF, just as in the case of quickly decreasing potentials, 
determines these functions in the unique way. 

The main result of the work is that, starting from a simple asymptotic characterization of EF in the weak sense, 
we get their asymptotic characterization in the traditional, uniform, pointwise sense. 

The possibility to use an explicit weak asymptotics for the description of the solution allows basically to determine a 
solution. On the other hand in the work [3] it is shown how one can use uniform asymptotic formulas for the numerical 
description of the EF. Besides, the result achieved offers quite a convincing hope for its possible development and 
derivation of the traditional complete proof that there exists a unique solution with the asymptotic behavior achieved 
here. It is not simple and in the present work we do not go beyond the heuristic construction of the asymptotic 
formulas for EF. 

Description of the model . Let us consider a more detailed description of the model. The initial configurational 
space of the system is R^. Having stopped the motion of the center of the masses, we arrive at the system on the 
configurational space F = {z : z e R^, z ~ {zi,Z2,Z3}, zi + Z2 + Z3 0}. On F there is a scalar product (z,z'), 
induced by the scalar product on R^. The system at F is described by the equation H'^/ = X'^, 5* ~ 'I'(z) e 
C, z S F, H = — Az + V{z), V{z) = w(xi) + w(x2) + ^(xa), Xj € R^. Here - is the Laplace operator at F, 

Xi = -^(z3— Z2),X2 = -^(zi— Z3),X3 = -^(z2— zi) . It is clcar that Xi +X2 +X3 = 0. Let us introduce also = y^Zj. 

It is easy to verify that at F yi +y2 +y3 = 0, z^ = (z,z) = (xj,Xj) + (yj,yj), j = 1,2,3 , A^ = A^ + Ay. 

Together with z e F, x, y G R'^ we will consider the dual variables, momenta q G F, k, p G R'^. We will assume 
that w(x) = 7, a > 0, though a generalization is possible for the case w(x) — ^ + w(x), xw{k) — >■ 0, x ^ 00. 

BBK - approximation. The plain wave at asymptotic description of the function ^(z, q) beyond the vicinities 
of the screens aj = {z e F, Xj = 0}, j = 1,2,3, must be replaced by the BBK-approximation \E'^^^(z, q). This 
approximation was studied in ^ , refer also to [6] , though it had been used also earlier in [71 18] . It acquires the form: 

vl/^^^(z,q) ^iVoe'<^'i>i?(xi,ki)i?(x2,k2)i^(x3,k3). (2) 

Here £'(x, k) — ^{—irj, l,ixk — i < x, k >), x, k S R^, 77 — $ - is a confluent hypergeometric function, refer 

3 , 

[3]. The constant A'o = H is a product of the normalization constants of the three two-body scattering states 

n!;'^ = (27r)-ie-^F(l + ir/j). 

It is noteworthy that the function 

Vc(x,k) =7V,e'<^'''>i?(x,k) (3) 

is the solution of the one particle scattering problem on the coulomb potential. 
The discrepancy BBK of the approximation Q[^^^^] can be easily calculated: 

g^^BBA'j ^ ^ k2-X2,k3-X3 > 'i'l'^^^-k^ki < ki-Xi,k3-X3 > <^[^2%-kik2 < k2-X2,kl-Xl > $'i$2*3, 

where ' defines the derivative by the last argument of the function $. With z — > 00 beyond some vicinities of the 
screens the discrepancy decreases quicker than the coulomb potential. Though the BBK-approximation is continuous 
up to the screens, its discrepancy on the screens decreases not quicker than the coulomb potential. Because of that 
the BBK-approximation cannot be used in the vicinities of the screens. 

Separation of variables. In the vicinity of each screen the Schredinger equation allows an essential simplification. 
For example, the complete potential V{z) = w(xi) + w(x2) +w(x3) in the vicinity of the screen cji is simplified because 
of the formulas X2 = — -^y — ^x, X3 = ^y — ^x, x = xi, y = yi- 

In the vicinity cti y >> 1, y >> x, therefore the formula V ~ v{x) + Vm{y), Vm — gives a good 

approximation to the potential. An equation with such a potential allows a separation of variables — A^x + {v{x) + 

Vm{y))x = Ax- 
Since we are interested in restricted solutions, for x naturally arises the following representation 

X(x,y,k,p)= y^,(x,k')^™(y,p')'5(fc''+y'-i?)i?(q,q')dk'dp', q-(k,p), q' = (k',p'). 
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Here ipc ~ is the solution of the scattering problem for the potential v, while ipm ~ is the solution of the scattering 
problem for the potential Vm- 

We can formulate now the goal of the present work more specifically. We are going to give an explicit description of 
the leading order of the asymptotic in the uniform topology relative to the angle variable, z. Precisely this problem 
is solved here on the heuristic level. Note that having this result, one can change the set-up of the problem (we assume 
that the vector q is beyond some small vicinities of the three screens) and start finding a solution, characterizing it 
by the uniform asymptotics. 

BBK— approximation in the vicinity of the screen. BBK-approximation in the vicinity, for example, of the 
screen txi but, however, in the domain where the discrepancy of this approximation decreases still quicker than the 
potential, naturally separates into a product 

^BSi^(z,q) =^c(xi,ki)^'i(z,q), *i(z,q) =7V^'')e'<''i^Pi>i?(x2,k2)i?(x3,k3), N^''^'> = N^:^'^ N^^l (4) 

In the vicinity of the screen the variables Xi, yi asymptotically have a different order, yi 3> xi. Let us compute 
the weak asymptotics in these assumptions ^i(z, q) for yi — ^ oo 



*i (z, q) = 5(p\ , ~y 1 )i3o (q) ^e-'^^P^ e'''^ [""^^ +^ " ""^1 e"'^' [^^^ 



(5) 



-'5(pi,yi)So(q) 
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Jyipi+iu Inyi iri2 In [z+ + i^V + j iria ln[z+ + ^|^y3+] 



Here we have used the following notations = ^(1 ± (pi,k2)). 



7± _ -/a 



(1 T (Pi,k3)), 



(xi,k2±pi), = (ii,k3 TPi), w = 2f7 + 2fT' 



i3o(q)=i?o(q) = -(2^)-^fcrfcr- 



Coefficients Z^^^y ^2(3) differ from the coefficients Z~^^y ^2(3) ^■^ change of the vector pi for — pi. 
Computation of the coefficient R. We will certainly need the asymptotic behavior x for y 
substitute the weak asymptotics tjjm into the integral x- 



00. Let us 



X 



V'c(x,k)— ((5(y,-p) 
py 



-ipy-\-irim In y 



Sm{y,p)e 



ipy — irjm In y 



) 6{k'^ + p'^ - E)R{q, q')dk'dp'. 



(6) 



Here 5*^ - is the coulomb scattering matrix corresponding to the potential Vm- 

For the further asymptotic simplification of the integral x for y — 00 we need the information on the structure 
of the coefficient R. It can be obtained from the comparison of the weak asymptotics on the variable y (the direct 
pointwise comparison turns out to be a significantly more complicated problem) of the integral (IH]) and of the BBK- 
approximation (|3][S]) in the domain, where the discrepancies of both approximations decrease quicker than the potential. 
It is a rather cumbersome computation that we are not able to demonstrate in the present paper. But it is the most 
crucial part of the work. The main idea of the computation is, however, quite natural. The representation of the 
BBK— approximation at some distance from the screen by the integral x is, in fact, the spectral resolution with 
the respect of the eigenfunctions ipm , and therefore, the resolution coefficient, i.e. the coefficient R, can be found 
explicitly. 



The result is: 



■s k 



, k + (p' - p)Bin 

\k+{p'-p)B 




(7) 



kpk'p' ip'-p + iO) 1+'^ 



\k+{p'~p)Bout\J 

Here the kernel G = satisfies the equation J 5m(y,pp')G'(p', p)(ip' = (5(y, p) 



2a 

V3p' 



-|jr(l-*a)e¥BS", 



2a a I a 

V3p ' ^ ~ 2k2 2ka ' 

The coefficient Ain and the vector Bi„ can be found from the equations Ain 

^ oA:^*^' (-^m, k) = 0. 

Correspondingly, the coefficient Aout and the vector Bout can be found from the equations Ao^t 

ib)e-^B^^\ Bout = fk-^nouu (B„„t,k)=0. 
We used here the notations 



-Bin 

-^r(i 



i?™(q) = (2^)- 



^(l-(P,k2)) 



IV2 



^(l + (P,k3)) 



ins 



1.1V2 IJ-Vi 

'>'3 , 



(q) = (27r) 



-3 



^(l + (P,k2)) 



")2 



^(1 

2 ^ 



-1 ir)3 



(P,k3)) 



"'2 "-3 , 
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1 / k2 -p , _k3+P_\ „ 1 / k2+p , kg 



vsyi-ipM "i + (p,k3);' i + (P,k2) 'i-(p,k3>; 

With such a choice of R the asymptotic in the weak sense formula © for x and the asymptotic formula (2])- ([5]) for 
'^f^^ for y oo coincide in the leading order on the intersection of the domains 

Voj = {z : Xj < yj, < 1/ < 1}, Vij = {z : xj > , 1/2 < ^ < v}. 

The next step is the derivation of the weak asymptotics of the integral x (O with y 3> 1 and arbitrary x, x <^ y 
and the reconstruction of the pointwise asymptotics by the weak asymptotics in the vicinity of the screens: 

X(z,q) ^ 7Voe'<"'''^i?(xi,ki)i?(x2,k2)i?(i3,k3). (8) 

This expression coincides with the expression © for m^^^ with the substitution 

1 ^ VS lVkV'c(x,k) 

73 1 . V3 , .lVk^c(x,k) 

"^ = -^-2" ^ "^ = — y + ^2 ^.(x,k) ' 

where the functions V'cCx, k) are defined in the equation Q. 

This result is the central one. Note that at large x the expressions X2(3) turn into the ordinary expressions X2(3), 
while the integral ([5]) turns into '^^^^^ Note as well, that beyond the vicinities of the pair directions of the 
forward scattering (xj , kj ) = 1 , j = 1 , 2 , 3 , our result corresponds to the result of the work [S] . The uniform in angles 
description of the leading order of the scattering problem solution asymptotics has been obtained for the first time. 
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Formulating of the result. The sets Vo r\ IJ Voj and Vi = IJ Voj cover F: F = Vb U Vi. Consider the 

J = i i=i 

separation of the unit 1 = (o{x,y) + Coj(a^:y)j subordinated to this covering. Let us assume that there, where 

the functions Coji Co differ from the constants, they depend on the ratio p = Inxj/lnyj. 
Let us form the expression 

Let us also define on F the function 

3 
i=i 

We believe in the correctness of the following main assertion: Function '^"^^ correctly describes the asymptotic 
behavior of the solution ^ in the leading order in the uniform topology. One can also verify that the discrepancy 
Q[v[/'"*] — — Az^I''^* + V'^'^'^ — A^"" decreases at z — oo quicker than the coulomb potential. 
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